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Abstract 

We perform a NLO analysis of polarized deep inelastic scattering data to test two different 
solutions to the so called spin crisis: one of them based on the axial gluon anomaly and consistent 
with the Bjorken sum rule and another one, where the defects in the spin sum rules and in the 
Gottfried sum rule are related. In this case a defect is also expected for the Bjorken sum rule. 
The first solution is slightly favoured by the SLAC E154 results, but both options seem to be 
consistent with the CERN SMC data. 



PACS numbers: 13.60.-r, 13.60.Hb, 13.88. +e, 14.20.Dh 



1 Introduction 



The EMC experiment on g\ (x) jj| was giving for the first moment of the spin dependent proton 
structure function g\{x) at <Q 2 >= WGeV 2 

r?= C g\{x) dx = 0.126 ±0.010 ±0.015, (1) 
Jo 

a value much smaller than the prediction of the Ellis-Jaffe sum rule || for the proton (we use 
a 3 = Ga/Gv = 1.2573 ± 0.0028 and a 8 = 3F - D = 0.579 ± 0.025 §), 

T? FI = — - — = 0.185 ±0.003. (2) 
i£j 2 18 w 

This has stimulated a considerable experimental activity in measuring the polarized structure func- 
tions (SF) g\{x), gf(x) {g{ He ) and gf(x) |, | g. 

From the theoretical point of view, the hypothesis has been formulated 0] that a relevant 
isoscalar contribution, related, in the chiral-invariant factorization scheme [||, to a large gluon 
polarization, contributes to Eq. (0) at 0(a s ), such as 



H(Q 2 ) = j- 2 (as + f + iaole 2 )) (l - 2^i) - ^AG(Q=) 



(3) 



This would explain the defect in the Ellis-Jaffe sum rule for the proton, without affecting the more 
theoretically founded Bjorken sum rule ||, 

QCD corrections play an important role in modifying the r.h.s. of Eq. (0) which is, strictly speaking, 
valid in the limit Q 2 — ► 0, and reads 

(r?-r?)(Q 2 ) = ^(i-^), (5) 

where only the 0(a s ) has been retained. The validity of the Bjorken sum rule is so universally 
accepted that, in the analysis of the experimental data on polarized SF [10, 11, 12, 13|, it is almost 
always assumed, rather than tested. There are few exceptions to this attitude, as for example the 
paper by Altarelli et al. [11|, where is evaluated to be 1.19 ± 0.09, and the l.h.s. of Eq. (||) at 
Q 2 = 4 GeF 2 is 0.177 ± 0.014, in agreement with the predicted value at 0(a s ), 0.189 ± 0.002. 

An alternative interpretation of the defect shown by Eq. (|l|) has been given jl4| by relating the 
defects in the Ellis-Jaffe sum rule for the proton and the Gottfried sum rule [|1| for the isovector 
unpolarized SF of the nucleons, 

i F 2 p (x)-F 2 "(x) 12 /Y - 1 

dx = - + - / [u(x) - d(x)\ dx = -, (6) 

X 3 3 Jo 3 
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to be compared with the measurement performed by the NMC experiment |lq] , namely 0.235±0.026. 
The defect in the Gottfried sum rule may be ascribed to a flavour asymmetry in the sea of the 
proton, which does not mean the breaking of the isospin symmetry since the proton is not an 



isoscalar. A flavour asymmetry of the parton sea in the proton was advocated many years ago [17], 

in connection with the Pauli principle and an essential asymmetry in the proton which contains 

two valence u and one valence d quark. The sign of the asymmetry is just the right one to account 

for the defect in the Gottfried sum rule, d — u > 0, while to reproduce the value given by NMC one 

should have 1 

d-u= [ [d(x) - u(x)] dx = 0.15 ± 0.04. (7) 
Jo 

The conservation of the SU(3) non-singlet weak currents gives 

Au val =2F, Ad val =F-D. (8) 

Since in the proton is the valence parton with the highest first moment (the quark model sum 
rules would imply vJ = 1 + F = 1.459 ± 0.006), it is reasonable to assume, inspired by the Pauli 
principle and Eq. ([/]), that it is just this parton which receives less contribution from the sea and 
to write 

Au vai = 2F + u - d. (9) 

This implies a defect of about 0.033 (0.008) in the Ellis-Jaffe sum rule for the proton (neutron) and, 
what would be very relevant, a defect of about 0.025 in the Bjorken sum rule. In this framework 
one should expect 

rp F D 2, - 

2 18 9 (1Q) 

F 2D 1 
r? = -- — + —(u-<f) = -0.033 ±0.004, 

in reasonable agreement with experiment. 

Our strategy, to make a test for these two different interpretations, will be to compare with 
experimental data the predictions of parton distributions, whose properties will be dictated by 
these two assumptions. Since the data extend on a wide Q 2 range, we need a method to solve 
the integro-differential evolution equations [18] at the next-to-leading order (NLO). To this aim, 
we choose to reconstruct the SF by using a truncated Jacobi polynomials expansion Jl^, 20]. The 
main advantage of this method is that it is fast, since one can analytically calculate the moments 



(Mellin transforms) of the SF which enter the expansion. Indeed, it was already used in Ref. [21] 



for a NLO fit of the unpolarized distributions and in Ref. [13] of the polarized ones 



The paper is organized as follows. In the next section we will present the parameterization used 
to describe the SLAC data, with the different options inspired by the two interpretations previously 
considered. Then we will present the Jacobi reconstruction method and recall the QCD evolution 
of the moments. Finally, we will report the results of our analysis for the SLAC data and the QCD 
evolution to the SMC higher Q 2 value, before giving our conclusions. 



3 



2 Description of the SLAC data 



With respect to our previous paper [22|, where the analysis of the data was performed at leading 
order in a s and in the chiral-invariant factorization scheme, here we take advantage of the knowledge 
of the NLO anomalous dimensions and coefficient functions within the MS renormalization scheme, 
which is the most usual one chosen in the gauge-invariant factorization scheme. Therefore, instead 
of Eq. (13) of our previous paper, we write (we assume a flavour symmetric sea, Au sea = Au sea = 
Ad sea — Ad sea — As sea — As sea = Aq sea ). 

xAu val (x,Ql) = rj u A u x a -»(l-x) K (l+j v x), 
xAd va i{x,Ql) = -q d A d x a v{l- x) bd {l + -i v x), 

(11) 

xAq sea (x,Ql) = xAq sea (x,Ql) = r} s A s x a °(l - x) bs (l + 7 s z), 
xAG(x,Q 2 ) = VG A G x a G(l-x) b °(l + lG x), 

at Qq = 4 GeT^ 2 , where rj q (q = u, d, s, G) are the first moments of the distributions and A q = 
A q {a q ,b g ,j q ) is given by 

- li d »"'- 1 < 1 - + = f 1 + ^^n) (12) 

in such a way that 

/ dxA q x a o- l {\- x) b o{l + ~j q x) = 1. (13) 
Jo 

We fix 7s = 7g = 0, since it is not possible to fit the values of these parameters with sufficient 
accuracy. 

We fit the asymmetries 

A« n \x,<?) = = 9 \\ {x > Q2) 2x(l + R^(x,Q 2 )), (14) 

which are just the quantities experimentally measured. The choice of A\ would also allow one to 
minimize the higher twist contribution (that should partly cancel in the ratio in Eq. fll4|)), instead of 
cutting the data at low Q 2 . In constructing the ratio in Eq. ( |l4|) , we use the MRS parameterization 
at Qq = AGeV 2 |23| for the unpolarized distributions. The deuteron asymmetry, Af(x,Q 2 ), is 
given by 

gf(x,Q 2 " 



with 



Ff-(x,Q 2 ) 

gf(x, Q 2 ) = \ (i - \^d) {g{{x, Q 2 ) + g?(x, Q 2 )), (16) 

where ujjj = 0.058 [24] takes into account the small D-wave component in the deuteron ground 
state. We only include in the fit the data obtained by the SLAC experiments (||. The HERMES 
data H are not used since their lower precision does not change the results. We will compare the 
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results with the SMC proton and deuteron data by evolving the distributions to the Q 2 of the 
CERN experiment. 

In all the cases we will fix 

Vd = F - D = -0.339 ± 0.013, (17) 

as expected from the conservation of the SU(3) non singlet weak currents. The two different 
interpretations mentioned above will be characterized, for the first option (see below, fits A, B, 
and C), by the additional constraint 

Vu = 2F = 0.918 ± 0.013, (18) 

in such a way to obey the Bjorken sum rule with the first order correction, Eq. (||), consistently 
with the NLO approximation. For the second option (see below, fits D, E, F, and G) rj u will be 
a free parameter and the Bjorken sum rule will departure from the theoretical expectation. The 
sea and gluon parameters, rj s and ijg, will be left free in both approaches (we impose rja < 3 since, 
for high values of this parameter, the \ 2 shows a weak dependence on it) and we will also consider 
cases with vanishing rj s or (and) rjc to evaluate from the data how much their contributions is 
needed in both options. We also consider a case (fit A) with r] s fixed by imposing that the NLO 
value of rf at <Q 2 >= 3GeV 2 should be the one reported by SLAC E143, namely 0.127 ± 0.011. 
Further limitations on the parameters are: 

a v > a s > ac 

(19) 

b u > 3.96, b d > 4.41, b s > 10.1, b G > 6.06, 
in order to satisfy the positivity constraints with respect to the unpolarized distributions. 



3 The Jacobi reconstruction method 



There are different methods to solve the evolution equations 18] for the distributions (or the SF). 



One possibility is to numerically solve them in the x-space [25, |26|]. A faster choice is based on the 
use of the Mellin moments, 

F n {Q 2 ) = [ x n ~ l F(x,Q 2 ) dx, (20) 
J o 

since they are analytically calculable for distributions like the ones in Eq. (|il"l). After transforming 
the integro-differential equations to the n-space, they become ordinary differential equations in the 
Q 2 variable, which are analytically solved, and the solution is then numerically inverted to give the 
distributions in the x-space. 

An alternative choice, which avoids the numerical inversion, is to express the SF by means of 



orthogonal polynomials [27]. On the one hand, this method is more favorable because it allows one 



to use the analytical expressions of the Mellin moments; on the other hand, a positive feature is 



5 



also their weak dependence on the values of the reconstructed function outside the region where 
data are collected. 

The expansion of the SF in terms of Jacobi polynomials, 0? (x), takes the following form: 



oo 

F(x,Q 2 ) = x a (l-xf £ a^\Q 2 )@f (x) (21) 

k=0 

where the Jacobi polynomials satisfy the orthogonality condition 

f 1 x a {l-xf Of ef dx = 6 kl , (22) 
Jo 

and the coefficients a^*^' '(Q 2 ) are defined by 

a ( f\Q 2 ) = f 1 F(x,Q 2 ) ef(x) dx. (23) 
Jo 

The expansion in Eq. (^Tj) becomes useful if we consider the power expansion of Jacobi polynomials, 

©f(*)=E #WW, (24) 

3=0 

with ((a)^ = r(a + k)/Y{a) is the Pochhammer symbol) 



c > ^ = ( j J i^+Ty, v na+k+Dnp+k+DM • (25) 



By inserting it in Eq. (031) , we get 



4 Q/3) (Q 2 ) = E cf\a,(3)F J+l (Q% (26) 
i=o 

and obtain for F(x, Q 2 ) 

oo k 

F(x,Q 2 ) = x a (l - xf Yl @ f( x ) E c f ^+i(Q 2 )- (27) 

fc=0 j=0 

This formula is very useful because the Q 2 dependence of F(x, Q 2 ) is contained in its moments, for 
which the solution of the evolution equations up to the NLO is well known. Evidently, one has to 
approximate Eq. ( p7[ ) by truncating the infinite series to a finite number of terms, 

N k 

fW(x) = x«(i - xf x: jfc*) E #W) ( 28 ) 

fc=0 j=0 

This truncation implies that outside the range where data have their values the approximation 
can fail, and it is possible that oscillations take place; however, by choosing a suitable value of N, 
one can reconstruct the function with sufficient precision. In our analysis, we reconstructed the 
unpolarized SF with N = 16, a = —0.99, and = 4.03, while for the polarized ones we used N = 8 
and allowed a and f3 to vary. The values found in the different cases considered below are reported 
in Table |l|. 
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4 Evolution of the QCD moments 



Given the Mellin moments of the charge-conjugation even (odd) non-singlet, singlet and gluon 
unpolarized distribution, at Qq, 

QlsniQl) = I x n - l Ql s (x,Q 2 ) dx, (r? = ±l) 
Jo 

S n (Q2) = f 1 x «-is( X) Qg) dx, (29) 



G n {Ql) = I x n - x G{x,Ql) dx, 
Jo 



their evolution with Q 2 is given by the renormalization group equations. In the MS scheme one 
gets, for example for the non-singlet distribution, 



<y (o 2 )-o" fc 2 ) n ° (Q2) i 2/30 



l + ^£ W i _a s WoJ / 7ivs 7ivs Pi 



4^ \ 2/?o 2$ 



2 



(30) 



The complete expressions for the evolution of the distributions can be found in the Appendix, 
and (the last one appears in Eq. ( |A. 1|) ) enter the unpolarized anomalous dimensions for the 

non-singlet and singlet operators at NLO, in the following perturbative expansion, 

= Os_ f ( )n (l)nr, 



7 a "" = g [^ >n + ^""J , (31) 

(r/ = ±1 and r/ = 1 for the non-singlet and the singlet respectively), whose expressions was calculated 
in |2^, |30|]. We also have 

a s NLO {Q 2 ) _ 1 fclnln^ 



47r /9bi"S ^Zn 2 ^' (32) 

/3 = ll-fn /! ft = 102 - f n/, 

which is the QCD coupling constant at NLO (the LO expression is obtained taking (3\ = in 
Eq. (|3^)). According to the number of active flavors n/, the value of A^/^ is modified. We use 
A-nlo = 0.2263, so to have a s (M§) = 0.118 [31]. From the matching of a s at the quark thresholds 
we get, for m b = 4.5 GeV, A^ ] LO = 0.326 and A§£ = 0.3797, for m c = 1.5 GeV. 

Analogous expressions to Eqs. (jH) and JSH) can be written for the polarized distributions, 
where the polarized anomalous dimensions can be found in [flOfj . The polarized non-singlet and 
singlet distributions are defined as follows: 

Qns = a 3 + y = -Au va i - -Ad va i, 

— = —-Au„„j + -Ad va i. ( 33 ) 



Qivs = -a 3 + — = --Au va i + -Ad 



3 3 3 



£ = a = Au vai + Ad va i + 6Aq, 



seat 
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Figure 1: The polarized parton distributions for fits A and D at Qq = 4GeV 2 . Note that, in the case of fit 
D, sea and gluons are vanishing. 



For the evolution of the polarized moments we used the fixed-flavour scheme [32, 1C] with nj = 3. 

Since the coefficients in Eq. (25) contain in the denominator higher and higher factorials as 
N gets large, we have to take into account the corrections to the anomalous dimensions described 
in Ref. [20|. They depend on the type of distribution we consider, even or odd (+ or — prescrip- 
tion, respectively), and consist in the substitutions in the anomalous dimensions reported in the 
Appendix. 



5 Data analysis 

By considering the SLAC data of the E142, E143, and E154 experiments Q, we get the results 
reported in Tables || and |3|. 

A tentative study of Tables || and |3| is very instructive. It shows that one is able to get sat- 
isfactory fits to the data with assumptions inspired from both the two different options we are 
considering, namely r\ u and % fixed to the values expected by SU(3) non singlet current conser- 
vation, with the Bjorken sum rule obeyed and r\ s given by the defect in the EJ sum rule (fit A), 
or only rjd fixed to that value and without sea and gluon contribution for the option inspired by 
the role of the Pauli principle (fit D). The first option is more successful in describing the neutron 
( 3 He) data (a total contribution to the % 2 of 7 instead of 20 for the 19 data points of E142 and 
E154) and deuteron data (a contribution to the x 2 °f 55 instead of 63 for 56 data points), while 
the second one is slightly preferable for the proton data (a contribution to the x 2 of 55 instead of 
74 for 63 data points). We also note that, with the first option, by taking rj s = we get a rather 
bad fit which is fit C. 
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Figure 2: The prediction of fits A and D with data are showed for proton, deuteron, and neutron asym- 
metries. The curves correspond to the fits at the Q 2 of the data points for proton and deuteron, and at 
Qq = 4 GeV 2 for neutron. The HERMES data Q are showed, but they are not used in the fits. 

Fits A and D differ in the x dependence of the valence parton distribution which is more 
singular at x — > (x~ ' 427 ) for fit A than for fit D (x~ ' 119 ). They also differ in the behaviour for 
x — ► 1, where one has a faster power decreasing of the d partons for fit D ((1 — a;) 6,50 ) than for fit 
A ((1 — a;) 4 66 ). It is appealing, from the point of view of the interpretation inspired by the Pauli 
principle, to find in the case of the fit D, Au va i = 0.741 ~ 0.918 — 0.150 = 0.768, as expected by 
the connection to the defect in the Gottfried sum rule, but it is worth to remark that the better 
description of the s He data by fit A seems more relevant than the better description of the proton 
by fit D. When we leave free r\ u and r] s , with or without gluons (fits F and G), we find slightly 
better fits, with an almost average verdict for r) u in fit F: 

rfi = 0.822 ~ !k±!k = 0.830, 

rjf = -0.040 c^-n*, ( 34 ) 
5 

Vg = L35 - Vg- 

In Fig. |] the polarized parton distributions are plotted for fits A and D and here we note that 
fit A is very similar to a solution obtained in Ref. [jl3[ (see their Fig. 5). In Fig. Q we compare the 
predictions of the same fits with the SLAC data. In Fig. || we show the evolution of fits A and D 
to <Q 2 >= 10 GeV 2 , compared with the SMC data, which are consistent with both of them. For 
this evolution we used the fortran code by D. FaschingQ p5| , because the SMC data extend on a 
different range of x and Q 2 , and the previously calculated a and (3 in the Jacobi reconstruction of 
the structure functions would not allow one to well describe them. 



*We thank the author for providing us his code. 
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Figure 3: The predictions of fits A and D for the proton and deuteron g\ structure functions at <Q 2 >= 
lOGeV 2 are compared with the most recent data obtained by SMC (see Ref. ||). 

6 Conclusions 

We have studied the most accurate available data on polarized deep inelastic scattering with the 
purpose of testing two possible theoretical interpretations of the spin crisis. As a result of our 
analysis, we can conclude that the interpretation related to the Pauli principle is not contradicted 
by the rather precise SLAC data which, however, slightly favour the generally accepted solution 
with the Bjorken sum rule obeyed and a sea contribution to the proton spin, in the gauge-invariant 
factorization scheme. Obviously, a better determination of the small x region could help to clarify 
further the present situation. 
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Appendix 



The complete solutions of the evolution equations for the non-singlet, singlet and gluon distri- 
butions, in the MS scheme, are: 



(Qo) 



a. 



(0)n 
. 1 NS 

NLO (Q 2 )\ 200 



OLs NL °{Ql)) 

£„(Q 2 ) = [(1 - On) S n (Q2) _ a n G n (Q 2 



1 | a s LO (Q 2 )-a s LO (Q 2 



47T 



, NLO {Q 2 )\ 2/3 ° 



v 2 A) 2/3 2 



«/ LO (^) 



47T 

a s LO (Q 2 ) 



+ 



X n -A™ 



(1)7 

7+- 



+ 



2/3o 2/3 2 
I- 2 ? • A': A" J ' 

A™ 

] f a s NLO (Q 2 )\ ^ 
a n Sn(Q ) - a n G n (Q ) 1 " 1 



4vr V«^ (^o) 



+ 



a/ (Q 2 )-a/°(Q 2 ) / 7^ A!^ 
4vr I 2/3 2/3 2 



+ 



A™ -A" 



a s LO {Ql) ( a s LO (Q 2 ) \ 2 "o 
4vr \a s L°{Ql)) 



+ 



47T 



(1)71 

2/3 + A™ - A" 



G„(Q 2 ) = [a n G n {Q 2 )-e n S„(Q 2 )" 



a, 



ATLO^Q2^\ 2/3 



a 



NLOtrfl 



(Ql), 



■ < 



1 a/ (Q 2 )-a/°(Qg) / A^ft 

4vr 2/3 2/3 2 ' " 



a s LO {Ql) ( a s LO (Q 2 ) 
4vr l^(Q 2 ) 



A™ -A^ 
2/3 



+ 



«/°(Q 2 ) 

47T 



(1)7. 

7+- 



2A + ^-A» 7 W«_An 



+ 



+ 



(1 - an) G n {Ql) - e n S„(Q 2 ; 



NLO (Q 2 \ \ 2 Po 



a s NLO (Q 



oc s NLO (Q 2 ^ 



0, 



(1)77 

7-1 



A!L/?i 



a/ (Q 2 )-a/°(Qg) 

4vr I 2/3 2/3 2 



+ 



4vr ^a s i0 (Q 2 ) 



A"-A n 



+ 



(A.l) 
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\Q 2 



in 



The quantities a n , a n , e n , AJ, 7±±, T^s and 7^," are denned and given in [33]. The corrections 
to the anomalous dimensions referred in Section |I| consist in the following substitutions [ 34 1 : 



(l)n 
7-+ 

2 /3 + A™ - A™ 



(i)n 



(0)n 
(0)n 



Ai 



(-1)" 
S(n) 



±1 



(-1)" {±^(§) + r? ± (n) [-25 2 (n) + C(2)]} , 
(-1)" {±^(f ) + r/±(n) [-45 3 (n) + 3C(3)] } , 



(A.2) 



-1) 



±S(n) + 77±(n)§C(3) 



where the series S'fc(n), S' k (n/2), and 5(n) are defined for example in Ref. and £ is the Riemann 
zeta function. 
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TABLE 





a 


P 


A 


0.173 


3.03 


B 


1.14 


4.20 


C 


0.265 


2.27 


D 


1.37 


4.32 


E 


1.34 


4.21 


F 


0.247 


3.45 


G 


1.36 


0.157 



Table 1: The values of the parameters a and /?, defined in Eq. (|2j), for the polarized structure functions, 
corresponding to the different fits (see section 3) are reported. 
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TABLE | 





A 


B 


C 


Vu 


0.918 ±0.013* 


0.918 ±0.013* 


0.918 ±0.013* 


Vd 


-0.339 ±0.013* 


-0.339 ±0.013* 


-0.339 ± 0.013* 




-0.065 ±0.019* 


-0.050 ±0.015 


0* 


m 


1.69 ±0.64 


0* 


3.00 ± 0.08 


ct v 


0.573 ±0.057 


0.587 ±0.057 


0.320 ± 0.046 


a s 


0.338 ± 0.080 


0.587 ±0.063 


- 


ac 


0.338 ± 0.082 


- 


0.312 ±0.106 


b u 


3.96 ±0.02 


3.96 ±0.02 


3.960 ± 0.008 


b d 


4.66 ±0.35 


4.71 ±0.35 


5.79 ±0.37 


b s 


10.1 ±0.8 


10.1 ±0.3 


- 


b G 


6.06 ±0.22 


- 


8.54 ± 1.08 


lv 


7.53 ± 2.06 


7.38 ±1.93 


28.3 ±7.7 


EJp 


0.128 ±0.012 


0.137 ±0.009 


0.167 ±0.003 


EJ n 


-0.061 ±0.012 


-0.052 ± 0.009 


-0.022 ± 0.003 


Bj 


0.189 ±0.003* 


0.189 ±0.003* 


0.189 ± 0.003* 


X 2 /NDF 


1.07 


1.07 


1.64 



Table 2: The values of the parameters found with the parameterization in Eq. ([Il]) for the different fits 
corresponding to the first option (see text for more information) are reported. The results for the r.h.s. of 
Ellis- Jaffe and Bjorken sum rules are also showed. The asterisk * indicates fixed values. 
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TABLE | 





D 


E 


F 


G 


'/U 


o 741 ± n ni6 


754 ± 021 


822 ± 024 


820 ± 038 


'Id 


_n qqo xnniQ* 

\J . *J (J %j _l_ VJ.U-L*J 


_n QQQ + f) 01 3* 




_n qqq + n ni 3* 


'Is 


0* 


0* 


—0 040 ± 020 


—0 028 ± 011 


Tin 

'IG 


0* 


3 00 ± 3 00 


1 35 ± 23 


0* 


n 


881 ± 087 


853 ± 091 


730 ± 075 


682 ± 090 


n 

Uig 






n von 4. n 1 fjo 


n 509 + n 220 






053 ±0152 


730 ±0185 




h 

u u 


tj • :J \J\) _1_ \J •\J\J\J 


3 Qfi 4- n 06 


3 96 ± 04 


3 960 ± 003 


u d 


6 50 ± 46 


6 25 ± 47 

\J , £i\J _1_ V/. J 1 


5 36 ± 44 


5 41 ± 50 

!_*.t:i _i_ yj.xjyj 


h 

u s 






41.5 ± 11.1 


10.1 ± 0.4 


b G 




6.06 ±3.91 


6.06 ± 2.48 




lv 


3.70 ± 1.56 


3.97 ± 1.66 


5.64 ± 1.89 


6.65 ±2.50 


EJp 


0.132 ± 0.003 


0.134 ±0.004 


0.124 ±0.013 


0.131 ±0.010 


EJ n 


-0.031 ± 0.003 


-0.030 ± 0.003 


-0.051 ±0.012 


-0.044 ±0.007 


Bj 


0.162 ± 0.003 


0.164 ±0.004 


0.175 ±0.004 


0.174 ±0.006 


X 2 /NDF 


1.05 


1.07 


0.97 


1.01 



Table 3: Same as in Table for the fits corresponding to the second option. 
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